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The high- Tc cuprates are possible candidates for d-wave superconductivity, with the Cooper 
pair wave function belonging to a non-trivial irreducible representation of the lattice point group. 
We argue that this d-wave symmetry is related to a special form of the fermionic kinetic energy 
and does not require any novel pairing mechanism. In this context, we present a detailed study 
of the bound states and resonances formed by two lattice fermions interacting via a non-retarded 
potential that is attractive for nearest neighbors but repulsive for other relative positions. In the 
' case of strong binding, a pair formed by fermions on adjacent lattice sites can have a small effective 

^\ I mass, thereby implying a high condensation temperature. For a weakly bound state, a pair with 

0^ ■ non-trivial symmetry tends to be smaller in size than an s-wave pair. These and other findings are 

' discussed in connection with the properties of high-Tc cuprate superconductors. 

H ■ 

, PACS numbers; 71.10.Fd, 74.20.-z, 03.65. Ge 



C/3 



I. INTRODUCTION 



^ ' Recently, through angle-resolved photoemission spectroscopy 0] and the tri-crystal Josephson effect ||] , a substantial 
: d-wave component has been detected in the order parameter of the cuprate superconductors. For tetragonal crystals, 
' such as Tl-based cuprates in which the square symmetry of the Cu02 plane is perfect, the order parameter appears to 
i: be entirely d-wave. In other superconducting cuprates, with orthorhombic distortions, the order parameter contains 
■ a large fraction of c?-wave. These observations have prompted theorists to look for exotic mechanisms that could lead 
\^ , to a non-trivial symmetry of the order parameter. It may, however, be true that it is simply the fermionic kinetic 
^\ ' energy that leads to d-wave superconductivity. An exactly solvable problem with two fermions would provide a strong 

, argument for this possibility. 
Q ■ Pines et al. have proposed a one-band model for cuprate superconductors in which quasiparticles hop between 
O I copper sites with amplitude ti < along the sides of square cells and with amplitude t2 > along the diagonals. 

I ■ The interaction between quasiparticles is mediated by magnon exchange and, in the static limit, can be thought of 
Oh as a repulsive potential for two quasiparticles on the same antiferromagnetic sublattice and as an attractive potential 
^ otherwise. In the present paper, we demonstrate that a positive t2 is actually a necessary condition for a two-fermion 
system to have a d-wave ground state. In cuprate compounds, a reversed (positive) sign of t2 is indeed expected if 
doped holes occupy dj.2_y2 orbitals of the copper ions. Such orbitals have four lobes along the principal axes of the 
CuO planes, with a positive amplitude along one axis and a negative amplitude along the other. For adjacent copper 
5^ ions the overlapping lobes have the same sign, while for next-nearest neighbors positive lobes overlap with negative 
ones. Therefore, a reversed sign of ^2 occurs. 

Measurements of the upper critical magnetic field have revealed a very short coherence length in high-Tc cuprates. 
As we shall see, for a d-wave pair, the average distance between the two fermions is on the order of the lattice spacing, 
even for weak binding. Therefore, these pairs are very localized and may well undergo a Bose-Einstein condensation. 
The resulting transition temperature Tbe decreases with a growing effective boson mass. For a point-like boson on 
a lattice (strong on-site binding) , the mass is large because the particle can only move via virtual fermionic states in 
which the pairing is broken As we shall show, this is not the case for an extended boson such as the d-wave pair 
discussed in this article, and the suppression of Tbe disappears. 

Phcnomcnological models with resonance bosons have been proposed by Friedberg and Lee |^ and by Micnas, 
Ranninger, and Robaszkicwicz The s-wave description given in j^] can be extended to include the possibility of 
resonances having non-trivial symmetry. In the present paper, we discuss the formation of a two-fermion resonance 
(s-wave or d-wave) and show that such a resonance boson can propagate long distances before decaying. 

The present article is organized as follows: In Section we describe the fermionic Hamiltonian and present a 
theorem concerning the symmetry of the two-fermion ground state. Deep bound states, shallow bound states, and 
resonances are analyzed in Sections III , [V , a nd |v|. The discussion in these sections is mostly restricted to the case of 



a square lattice, with the exception of Section III A , in which we deal with a face-centered cubic lattice as in fuUerene 
superconductors. Concluding remarks are made in Section 
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II. TWO FERMIONS ON A PERIODIC LATTICE 



A. The Hamiltonian and its properties 

We consider a periodic lattice of M sites with symmetry under inversion at each site. The ferniionic kinetic energy 
is represented by a sum of "hopping" terms, 

Ho=Y.t{v-v')al{v')a,{v), (2.1) 

s,r,r' 

with s = t (spin up) and | (spin down). The interaction between two particles is described in terms of a potential 
energy that is a function of the relative separation and the total spin of a pair: 

i/i=^[g5(r-r') 5(r, r')^5(r, r') 

r,r' 

+ 5T(r-r') ^r;(r,r')r^(r,r')], (2.2) 
/i 

where (7s(r) and gxi'^) represent interaction energies in singlet and triplet states with a relative separation r and 

where 5(r, r') and 7^(r, r') are the destruction operators for spin-singlet and spin-triplet states on sites r and r', that 
is, 

5(r, r') ^ «T(r)«i(rO + «T(rOai(r) ^ (2.3a) 

To(r, r') = «T(r)«i(r')^«T(rOai(r)^ ^^.Sb) 
^ , t-. at(r)at(r') 

r+(r,r') = ^^^^i^, (2.3c) 

T_ir,r') = ^^^^}^. (2.3d) 
V 2 

The particular normalization of these operators has been chosen for later convenience. In this paper we are mainly 
concerned with singlet states. Therefore, the triplet part of the interaction will not contribute, and the subscript in 
the interaction energy gs will usually be suppressed. 

It follows from ( ^.l| ) that the one-fermion energy spectrum is 

ep = ^t(r)cos(p-r) (2.4) 

r 

with p the Bloch momentum. The corresponding effective mass in D dimensions is given by 

l=-lj2tiry (2.5) 

r 

for the case in which the inverse mass-matrix {■m~^)ij — 6ij/m. In this paper, we only consider the case for which 
m > 0. 

A convenient basis set in the Hilbert space of two fermions consists of states with total Bloch momentum P, relative 
separation r, and spin projections si and 

|P,r,si,.2) = -^^e^P-^-^'+r/^) aUr')al{r' + r)\0). (2.6) 

Applying the kinetic energy operator ( ^ ) to state (2^), we obtain 

Ho\P, r, si, S2) = (P-i-'/2)|P, r + r', si, 52)- (2.7) 
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The problem of two fermions is therefore reduced to that of a single particle hopping over r sites (in the relative 
position space) with an amplitude 2t(r) cos (P-r/2). Eigenvectors of Hq are states with total momentum P and 
definite relative momentum p, 

|P, p, si, S2) = a(P/2)-p,.i a(P/2)+p,.2 (2-8) 

with eigenvalues 

Ep.p = e(p/2)+p + e(P/2)-p 

= 2^t(r)cos(Pr/2)cos(p-r). (2.9) 

r 

The operators ap^s and as{v) are related by a Fourier transformation, 

l=^e-*P--a,(r). (2.10) 



The projection of any two-fermion state on the state ( ^.6D is the wave function in relative position space: 

V'(r) = (P,r,si,S2|P,V',si,S2). (2.11) 

Plane waves, 

^ ^(r|p) = (P,r,si,52|P,p,si,52) = ^e^P■^ (2.12) 



form a complete orthonormal set of functions: 

l5]V'*(r'|p)V'(r|p) = 5r'r, (2.13) 



p 

-i^V(r|p')V'*(r|p)- Vp- (2.14) 

r 

For singlet states, the two-fermion wave function must be symmetric under the transformation r — r: 

^ V^(r|p) = (P, r|P, p) = J|p cos (p-r) (2.15) 



for r ^ 0. For two fermions on the same site (r — 0), '0(O|p) = 1 is already symmetric. Using the symmetrized wave 
functions ( ^.15 ), we can write the interaction matrix elements of Hi as 



(P,q|i/i|P,p) = -i^V*(r|q)5(r)^(r|p). (2.16) 

r 

Sites r and — r are now considered equivalent and only one of them enters the sum. 

B. Schrodinger equation 

To solve the Schrodinger equation, 

{Ho + Hi)\P,^) =E\P,^), (2.17) 
for a two-particle bound state, we write |P, ■0) as a sum over states having definite relative momenta: 

|P,^) =^F(p)|P,p). (2.18) 
p 

The Schrodinger equation is then transformed into an integral equation for F{p): 
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{E - £;p,p)F(p) = Pl^i |P' (2.19) 

q 

= -l5]V/(r|p).g(r)^F(q)^(r|q). 

r q 

Recalling that V'(r) — Sq ^(fl)V'(r|q), we obtain the foUowing equation for the wave function: 

V'(r) = ^(^' rOsCrOV^r'), (2.20) 
where the Green's function G is defined by 

q 

In what foUows, we shaU truncate the interaction potential so that g(y) vanishes beyond a certain distance, that is, 
g{v) 7^ only when r = 0, ±Ri, . . . , ±R„„i. Equation ( 2.2C| ) for these sites can then be cast into matrix form: 



[l-g{E,V)g]i^ = Q, (2.22) 

where ip is an nxl matrix with t/'j = ''I'O^j), QiE^V) is an nxn matrix with CJ.y(i?,P) = G (E, P; R^, R,-), and g is 
a diagonal nxn matrix with ga = gCYLi). Note that the kinetic energy of two particles ( |2.9| ) can also be written in 
matrix form, 

£^P,P = e(P/2)+p + e(P/2)-p 

= 2^V*(P/2)r,,V^,(p), (2.23) 



where the diagonal matrix r is defined by th = i(Ri) and where V'j(p) = ■0(R-i|p)- It follows from Equation (2.22) 
that the eigenvalues E are determined by the secular equation, 

det [^;(£;,P) = 0. (2.24) 



For each eigenvalue E, the wave function at th e n sites V'(R-j) is determined by Equation ( 2.22 ), while for sites beyond 



the interaction range one uses Equation (2.20) to determine '0(r) 



In the case of a scattering process, the solution to the Schrodinger equation (2.17) with an incident plane wave 
'(/'(rip) satisfies 



^(r) = V'(r|p) + Y GiE, P; r, r')5(r')V'(r')- (2.25) 



ip{r) is therefore given by 



V;(r)=^(r|p) (2.26) 

n n 

+ E E ^(^' ^' ^i)Aj{E, P)i;{R, |p), 

where the nxn scattering matrix A{E, P) is defined by the equation^ 

A-\E,P)=g-^~g{E,F). (2.27) 

A{E, P), as a fmiction of complex E, has a cut along the real axis corresponding to the continuum spectrum of two 
free fermions. Depending on the potential g{r), it may also have poles on the real axis (bound states) or below the 
upper edge of the cut (resonances). 



^It should be noted that Aij 7^ ^(r, r') |r=R; ,r'=Rj , where A is the scattering matrix defined on the entire lattice. It is true, 
however, that A~^ij = j4~^(r,r')|i.=R.,r'=Rj- 
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C. Symmetry of the ground state 



We shall now demonstrate that the wave function of the ground state of two fermions does not change sign unless 
some of the hopping amplitudes t(r) are positive. The situation is similar to that in the continuum, for which the 



ground-state wave function has no nodes. Note that a positive effective mass [Equation (2^)] requires at least some 
of the hopping amplitudes to be negative. 

We have seen that the problem of two fermions reduces to a one-body problem with modified hopping amplitudes. 
It therefore suffices to study the ground state of a single particle moving on the same lattice in the external potential 
g(j). We consider the case of total momentum P = 0, for which the amplitudes double (reminiscent of the reduced 
mass effect in the continuum). The ground state is obtained by minimizing the expectation value of the energy: 

(£;) =^2t(r-r')V*(r'Mr)+5]g(r)|V(r)p. (2.28) 

r,r' r 

Suppose th at t{ r) < for all r and that ip{r) is positive on some site ri and negative on another site r2. Then the 



first term in ( 2.28 ) can be lowered if instead of a negative V'(r2) we use |V'(r2)|, the second term remaining unaffected 
by such a change. Therefore, the absolute minimum of (E) must have V'(r) > on all sites, provided there always 
exists a sequence of allowed hops from any one site to any other site on the lattice. 



III. DEEP BOUND STATES 



When the magnitude of an attractive potential exceeds the kinetic energy of two fermions, compact pairs are formed 
and can be treated as small-size bosons. Pairs formed by two fermions on the same site, however, have a vanishing 
mobility in the limit of strong binding. This strongly decreases the condensation temperature ^] , in contrast to the 
continuum case. It is desirable to have a model which does not have this lattice artifact. The pairing of fermions on 
adjacent sites is one way around the difficulty. We address this issue by working out two specific examples of lattice 
fermions having a strong nearest-neighbor attraction: first on a face-centered cubic lattice and then on a square 
lattice. 

The advantage of extended-boson models is offset by a new problem, which is also a lattice effect. While only 
two fermions can be involved in on-site attractive interactions due to the Pauli principle, extended-range attractive 
potentials can result in the formation of many-body bound states. However, quantum statistics can help avoid this 
type of phase separation because the interaction strength in symmetric (spin-singlet) and antisymmetric (spin-triplet) 
spatial states can be quite different. Let us assume that the interaction is attractive in a singlet state and repulsive in 
a triplet. Because a fermion cannot form more than one singlet bond at a time, it will attract only one other fermion. 
More exactly, quantum spin fluctuations will produce a fluctuating potential between fermions. Only in pairs are 
these fluctuations completely suppressed, and a state with pairs may be favored over other many-body states. For 
details, see Appendix ^ 



A. Deep bound states on a face-centered cubic lattice 

The following simple model describes tightly bound fermion pairs with a small effective mass (close to that of 
individual fermions). The small mass results from fermionic pairs which move through the lattice without having to 
break the bond. Triangular and f.c.c. lattices offer such a possibility. 



The Hamiltonian [Equations (2.1) and ( p.2| )] for the f.c.c. model contains nearest-neighbor hopping with amplitude 
i, on-site repulsion go^ and attraction gi for fermions on adjacent sites. On an f.c.c. lattice, there are twelve nearest 
neighbors r' around any one site r, with 

r'-r = ^(±x±y), |(±y±z), ^(±z±x). (3.1) 

X, y, and z are the unit vectors along each side of the cube and a is the length of each side. 

If Iffil ^ 1^1 J so that the bond between fermions is strong, the relative-position wave function 'il'{r) ^ only on 
the 12 sites surrounding the origin. We therefore expect 12 bound states with the same potential energy gi and, in 
general, unequal kinetic energies at a total Bloch momentum P 7^ 0. The kinetic energy is determined by the hopping 



amplitudes, 2< cos (P-r/2), as discussed in Section II A. 
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At P = 0, all hopping amplitudes are 2t and the energy levels have degeneracy (see Appendix there are a singlet 
Ai, a doublet E, and a triplet T2 (Figure p. As long as i < 0, the ground state for P = is Ai, with ^p{r) = l/-\/T2 
on the 12 sites around the origin and zero otherwise. From each of these 12 sites a fermion can hop to 4 of its nearest 
neighbors (see Figure 0(a)). For P 7^ 0, the boson dispersion relation to lowest order in P is given by expanding 
cos(P-r/2) to second order: 



E{P) =gi + 8t 



12 



E 



P-(r-r") 



(3.2) 



where r" denotes the four nearest-neighbor sites shown in Figure |l|(a) . Because E is an isotropic function of P to this 
order, we average over all directions of the total momentum and substitute |r — r"| — a/\/2 to obtain 



£;(P) =gi+8t 



tP^a 



2„2 



(3.3) 



which gives the effective mass of a boson as 



■nib 



3r 

ta2 



6mf. 



(3.4) 



B. Deep bound states on a square lattice. 



On a square lattice, the mass of a strongly bound pair will be finite if one takes into account next-nearest-neighbor 
fermion hopping (with amplitude in addition to nearest-neighbor hopping (with amplitude ti). Suppose the 
interaction potential has strength go > for particles on the same site, gi < for nearest neighbors, 52 > for next- 
nearest neighbors, and otherwise. For spin singlets, we have five distinct relative positions with non-zero interactions: 
Rq = (0,0), Ri = (0, ±a), R2 = (±a, 0), R3 = (±a,±a), and R4 = {±a,^a). In this basis, t = diag(0, ii, fi, i2, i2) 
and V — diag((7o, (71, (72j 52)- According to Equation (2.5), the fermion mass is 



nif 



2(-ti - 2t2)a2' 



(3.5) 



As in the determination of a single fermion mass, we expand the energy of a two-fermion bound state in powers of 
P. Near P = 0, it is convenient to use the irreducible representations of the symmetry group (see Appendix ^); and 
we choose the following five components ipi{p), 



V'(P) 



/Mp) \ 
Mp) 

-02 (P) 
V'3(P) 
V^4(P) / 



( 



1 



COS Pa; Q + cos p J, a 

COS Pa; a — COS p J, a 
2 COS Pa; a COS Py a 
y 2sinpa;a sin p^a J 



(3.6) 



as a basis set. Functions ipOj V'li s-^d "03 belong to the Ai representation (refered to as s-wave in the literature); ^2 
transforms as Bi (da;2_y2-wave); and ■04 as B2 (da;y-wave). Assuming that the nearest-neighbor attraction gi is strong 
and E gi, we expand the Green's function 



V'»(P)V'-(P) 



Gy(£;,p) = ^y — , 



(3.7) 



in powers of 1/E and thus determine the dependence E{P). Because 



(3.8) 



we need only the diagonal matrix elements to calculate the lowest-order correction to E. We have 
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E 



O 



Gn(£;,P) = ^ + ||^3(P/2) + o(ij 



i_||,3(P/2).0(J, 



Upon substituting these expressions into equation ( ^.24 ), we obtain the foUowing dispersion laws: 

1 



E = gi±2t2M^l2) + 



91 



(3.9) 



where the plus sign corresponds to the Ai state V'l and the minus sign corresponds to the Bi state ?/'2. The other 
three bands have no dispersion to this order. For a positive (negative) ^2, the ground state transforms according to 
representation Bi {Ai) and the mass of the boson is determined strictly by the diagonal hopping term: 



-2<i - Uo 



rub 



\t2\a^ 



\t2\ 



-ruf. 



(3.10) 



Note that, unlike in the case of a strong on-site attraction, an extended boson can easily move through the lattice 
without breaking its bond, hence a small mass even for a strong attraction. It can even be smaller than the mass of 
its constituent fermions. 

It is straightforward to see what determines the symmetry of the ground state. The relative position of two fermions 
is restricted to four nearest-neighbor sites by a strong attraction. Depending on the sign of t2 , the wavefunction must 
either alternate its sign on these sites or maintain its sign to minimize the kinetic energy. 



IV. SHALLOW BOUND STATES 



In this section, we investigate the properties of bound states that lie close to the bottom of the two-fermion band, 
which is assumed to have a single minimum Eq at the relative momentum p = 0. As described in Section II B, we first 
solve Equation (2.22) on the relative sites having a non-zero interaction potential. As before, we assume that g{r) ^ 
only for a finite number of lattice sites r = 0, ±Ri, . . . , ±R„_i. Let us also assume that the same-site repulsion is 
infinite so that "0(0) = 0. Using Greek indices to enumerate sites ±Ri, . . . ,±R„_i, we obtain from Equation ( 2.22| ) 
the following set of equations for the vector A = 51/) :0 



^OoAo + Goa^a — 9 ^QQ^Q — 0, 
GaO^Q + Gafi^p — g~^a0^P- 

Summation over repeated indices is assumed. Eliminating Aq, we obtain 

GaoGof) 



^ GaoGaf}\ A -1 A 

Gall A/3 = 5 afjA/^ 

yoo / 



Now we can use Equation (2.20) to calculate the wave function on the remainder of the lattice: 

^(r) = G(r,R,)A, 



G(r,R„) 



G(r,0)G(0,R„) 



G(0,0) 



A, 



(4.1a) 
(4.1b) 



(4.2) 



(4.3) 



^The BCS gap A(r) is the second-quantized counterpart of this function. 



7 



Apart from a multiplicative constant, the long-range behavior of tp{r) is determined by the properties of the Green's 
function at large separations, which in turn is controlled by the long-wavelength behavior of the summand in Equation 



(2.21). The latter diverges at q = as energy E approaches the bottom of the two-fermion band. Then the 



denominator in (2.21) can be approximated using the effective fermion mass: 



E - iJp., ^E-E^,o-^^ ^l+m^ (4.4) 

rrif ruf 

which defines a coherence length at themperature T = as the size of the bound state. For r 3> Co, the propagator 
decays exponentially with the distance, while in the intermediate range, a <^ r <^ ^q, it decreases as a power of |r — r'|. 

It is possible to isolate the divergent part of the Green's function in the form of G{E, P; 0, 0) = G(0, 0). We define 
a regularized Green's function as follows. For r 7^ and r' ^ 0, 

F(r - ^ V ['^'^'"''^^ ~ V'(r|0)][^*(r1q) - V^*(rlO)] 

f(..o)^iE '''"'"";!?"'"'""" . 



^ V(0|q)[V^(r'|q)-V^*(r'|0)] 



q 



These sums are regular at q = 0. By inspection, 



G(r, 0) = V(i-|0)G(0, 0) + F(r, 0), (4.6a) 

G(0,r') = G(0,0)V'*(r'|0)-hF(0,r'), (4.6b) 
G(r, r') = V(r|0)G(0, 0)V'*(r'|0) + V^(r|0)F(0, r') 

+ i^(r,0)Vj*(r'|0) -|-i^(r,r'). (4.6c) 



Upon substituting these functions into Equation ( |4.3D , we finally obtain the wave function as a sum of two terms, 

V<(r) = F(E,P;r.R„)A<, 

-mP;0,R„)A,||||^. (4.7) 



Consider states with total Bloch momentum P = 0. The first term in Equation (4.7) is proportional to the regularized 
propagator F{E, P; r, Ha), whose Fourier transform remains finite as q 0. In the intermediate range, a <C r <C Co, 
it decays with r as or faster (D is the number of spatial dimensions). The second term, on the other hand, is 
on the order of r^~^, as can be inferred from its singular Fourier transform. It therefore should dominate at large 
distances. However, this term vanishes if the state belongs to any irreducible representation other than the trivial one 
due to the prefactor F{E, 0; 0, Ra)AQ.0 This means that a bound state which transforms according to a non-trivial 
irreducible representation of the point group will be smaller in size than a symmetric bound state of the same energy. 
This may have verifiable consequences for superconductivity. 

To illustrate the importance of this point, we use a continuum approximation to the two-particle problem. In 
the intermediate region, a<Cr<?;Co, the two-particle wavcfunction 4>{r) satisfies the continuum Laplace equation, 
V^'(/'(r) = 0. In 2 and 3 dimensions, solutions with angular momentum / and its projection m are 

^(r) = e^-^r"™ (2 dimensions), (4.8a) 
^(r) = Yirn{9,4>)r-^-'^ (3 dimensions). (4.8b) 



^The Green's function F{E,0;0,Ila) is invariant under all point group transformations. 
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It is generally thought that the lattice spacing a is irrelevant as long as we deal with a weakly bound pair (^o ^ a). A 
simple scaling argument demonstrates that this is not always the case. To imitate the infrared and ultraviolet cutoffs, 
we set "0(1") =0 outside the interval a < r < £^q. In 3 dimensions, the n-th moment of r is then given by 



J a 



^0 2/-2+n^2^ 



(4.9) 



For / = one obtains 



(^"> = (4.10) 

for any n > 0. Indeed, an s state is not sensitive to the ultraviolet cutoff (the lattice spacing). For / > 1, 

2/ — 1 

n<2/-l: (r")^—- a", (4.11a) 

2« — 1 — n 

71 = 2^-1: (r") = (2Z- 1) a^'-^log^, (4.11b) 



71 > 2^ - 1: (r ") = 



21-1 / a 



n-2; + l ^" VC( 



a 

21-1 



This suggests that, for a shallow bound state with non-zero angular momentum, the measured size of a pair is 
determined by both length scales, ^^^d therefore may be much smaller than ^q. The physical reason behind 

this difference is that, at the same value of kinetic energy, states with a larger angular momentum have a smaller 
radial momentum. Therefore, the ability to tunnel into the classically forbidden region (to distances on the order of 
^o) is substantially reduced. Similar arguments hold for the two-dimensional case. 

The small size of bosonic pairs may imply a high upper critical magnetic field. As has been argued in |Q, the upper 
critical magnetic field for the condensate of an ideal charged Bose gas at zero temperature is infinite, because of the 
very different roles played by the Coulomb interaction in Bose and Fermi systems. Therefore, a rapid increase of 
the critical field as T — > should be expected for an ideal charged Bose gas. Recently, a similar behavior has been 
observed in Tc-suppressed cuprates 1^ . This observation indicates that Bose condensation may be dominant in these 
superconductors. Of course, bosons in real systems are not point particles. A rough estimate of the upper critical 
field would be the magnetic field at which the cyclotron radius of fermions is comparable to the size of d-wave pairs, 

(/)o - 2^i?c2(r'), (4.12) 

where 00 = hc/2e is a flux quantum. Formally, this criterion reproduces the Landau-Ginzburg result, provided that 
the pair size is replaced by the coherence length ^. We stress, however, that in a d-wave superconductor the former 
can be smaller than the latter. 



We now investigate the model described in Section [II B in order to determine the symmetry of a shallow ground 



state. Consider the case in which the bound state first appears below the continuum band. Using the basis set (p.q), 
we write the secular equation in the form of Equation ( |4.2| ). It is assumed that the on-site repulsion go is infinite. 
For a state with an energy just at the bottom of the two-fermion band, it is convenient to use a regularized Green's 
function T{E, P) as was introduced previously. Equation (|4.2|) then reads 



T^p-^^]^p^g-\,^p. (4.13) 



00 



t/oo becomes infinite as E ^ 2eo, which further simplifies the problem. The critical values of the coupling strength 
arc then given by the equation 

- 9-\p) = 0, (4.14) 

where Tap is to be evaluated at P = and E — 2eo. For the irreducible representations Ai, Bi, and B2; one obtains, 
respectively, 

(^11 - 5r')(-^33 -92')^ -^13^31, (4.15a) 

^22 -9i'^ 0, (4.15b) 

jF44 -£,-1 = 0. (4.15c) 



9 



As long as 32 represents repulsion, the B2 bound state (dxy) does not exist. 

If there is no next-nearest-neighbor repulsion (52 = 0), it can be seen that the ground state is Ai (s wave). Indeed, 



Equation ( |4.15a| ) reduces to Tn — and one obtains the following relation between the critical values of gi for s 
and dx2-y2 waves: 

Thus, it takes a weaker attraction to produce an s-wave ground state when there is no next-nearest-neighbor repulsion. 
When the next-nearest-neighbor amphtude ^2 is negative, the ground state is an s wave, in accordance with the 



result of Section II C . Therefore, the search for a non-trivial symmetry of the ground state is restricted to the domain 
in which both ^2 > and 52 > 0. The results can be summarized as follows. An increasing next-nearest-neighbor 
repulsion 32 raises the kinetic energy of an s wave but does not affect that of a dj.2_y2 wave. A positive i2 means 
positive kinetic energy for s waves and negative kinetic energy for d waves. The actual "phase diagram" is presented 
in Figure ||. For the value of t2 — 0.45|ti| inferred by Pines from fitting NMR results we find that the strength 
of the next-nearest-neighbor repulsion must be at least 0.35|ii| in order for the ground state to have the (x^ — j/^) 
symmetry. 

V. EXTENDED RESONANCE STATES 
A. General approach 

Even if the two-particle potential energy is positive at all separations, there still remains the possibility of two 
fermions forming a resonance state - a state bound for a limited time - provided that the potential has a local 
minimum. It has been shown that resonance bosons can form a condensate, thus inducing superconductivity in 
the boson-fermion system. 

On a lattice, we may approach the problem of a single resonance in the following way. When a local minimum of 
the potential is surrounded by a closed wall of sites with infinite repulsion, the wave function cannot leak out and two 
fermions form a stable bound state. Then the Schrodinger equation, {Hq + _ffi)|P, tp) — E\P, has a solution for a 
real E inside the two-fermion energy band. The wave function ip{r) is non-zero only for separ ation s r within the wall. 



which makes the Schrodinger equation readily solvable. Its solution also satisfies Equation ( 2.2C ), which determines 
the location of the pole of the scattering amplitude. In the limit of a stable state, the residue of the pole vanishes, 
reflecting the fact that incident free particles cannot penetrate the barrier of infinite repulsion. Lowering the height 
of the wall to a finite value makes the localized two-fermion state unstable, thereby shifting the pole away from the 
real axis into the unphysical sheet. In what follows, we determine the mass and resonance width of a bosonic state to 
lowest orders in the inverse height of the wall. 

Let the wave function of the stable state be ip{r), normalized to 1. Outside the wall, ip{r) ~ 0. ip{r) satisfies 



Equation ( 2.2C ), which now must be modified due to the existence of a branch cut of the function G(£',P;r',r) on 



the real axis of E: 

^l;{r) = G{E ± iO+, P; r, r')g(r')V'(r') (5.1a) 

r' 

Because ip{r) does not contain incoming or outgoing free particles, it is invariant under time reversal and therefore 



satisfies this equation on both edges of the cut. The complex conjugate of ( 5.1a ) is 



rir)=J2rir')gir')GiETM+,P;r',r). (5.1b) 

r' 

As before, let us assume that the inter-particle potential ^(r) 7^ only at the 2n — 1 sites r = 0, ±Ri, . . . , ±R„_i. 
We introduce a new wave function A(r) = g{r)ij}{r), which remains finite at a site r even as g{r) ^ 00 as long as the 
state has finite energy. Using the matrix notations of Section [IB, we rewrite Equations 5.1a and |5.1b in terms of the 
nx 1 column vector A: 

[g-^ -g{E±iO+,P)]A = 0, (5.2a) 
At [5-i-g(^T«0+,P)] =0. (5.2b) 
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In what follows, we will omit the variables E and P unless their presence is necessary to avoid ambiguity. 

By regarding the diagonal matrix as an independent variable, we can treat the position of the pole i5 as a 
function of g~^. When g~^{r) = on the walls, E is readily determined from the Schrodinger equation for the interior 
of the wall. By changing g~^ infinitesimally, we can trace the trajectory of E in the complex plane. Differentiate 
Equation (5.2a): 

(5-1 - g)dA + fdg-' - ^dE^ A = 0. (5.3) 



Multiplication of Equation (5^) by A^^ eliminates the first term because of Equation ( 5.2b ). Then, 

A^^AdE = AUg-^A. (5.4) 
dE ^ ^ ^ 

The factor in front of dE is merely —1 (see Appendix 

We thus have the following expression for the first-order displacement of the pole: 

dE = ~A^dg-^A, (5.5) 

which is real. Formally, dE can be written as ip^dgijj, the first-order term in the standard perturbation series, although 
the latter is not well defined for an infinite repulsion. 

Next, from Equation ( pT^ ) we obtain the first-order correction to A: 

dA - -A (^dg-^ - ^dE^ A 

= -A{dg-^g + gdE)il;, (5.6) 

where A~^ = g~^ — . Here, again, the substitution of dQ /dE by —Q^ requires justification (Appendix 

To obtain the second-order correction to the energy of the pole E, we differentiate Equation ( |5.3| ) once again and 
multiply by A^/2 to eliminate the term with d^A: 



lAt f + '^d^E\ A 

2 \dE^ dE 



+ At ( dg-i - ) dA = 0. 



(5.7) 



The second-order differential of the independent variable g ^ is zero. We simplify this expression by using 
A'^{dg/dE)A = -1 and ^ A'^{d'^g/dE'^)A = 21/'^^;^ and by substituting dA from (fj): 

y^E = ^^g^p + (g dg-^ + gdE)A {dg-'g + g dE)^p 

^ -tp^g-^^ + d(p^Ad(p, (5.8) 

where 

d(p EE {dg^^g + g~^dE)i! = {dg~^ - AUg-^Ag-^)A. (5.9) 



In Equation ( |5.8| ) we used the definition A^^ =g^^ — g and the fact that dg'^ip — (the potential has been changed 
on infinite-repulsion sites where the wave function vanished). The first term in Equation ( |5.§| ) is real. The imaginary 
part of the energy is thus given by the second term: 

hmd^E ^lmd(f>'iA{E±iO+,P)d<l), (5.10) 



*This relation can also be proved along the lines of Appendix ^ For higher derivatives, however, one will have to use the 
values of the Green's function on the entire lattice. Therefore, this type of analysis only works through second order. 
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We now derive an expression for the mass of a resonance state. The problem of a resonance formed by two fermions 
on the same site has been addressed in an earher work by Friedberg, Lee, and Ren They found that, as the 
Ufetime of such a resonance increases, its mass grows without Hmit. This can be understood as follows. Two fermions 
constituting a long-lived resonance must be located on the same site. However, the motion of such a boson is mediated 
by consecutive hopping of the fermions to an adjacent site. This means that the fermions must part and then reunite: 
the boson decays into two fermions and then reappears at an adjacent site. As the decay amplitude vanishes, so 
does the kinetic energy of the resonance, thus giving an infinite mass. The same problem exists in the case of bound 
states The situation changes if we consider a resonance formed by fermions on two adjacent sites. Now a pair of 
particles can move without changing the distance between them. On a triangular or face-centered cubic lattice, this 
can be accomplished with only nearest-neighbor hopping. On a square or simple cubic lattice, fermions must be able 
to hop to next-nearest-neighbor sites. Under such circumstances, even long-lived resonances will have a finite mass, 
which can readily be determined. In the limit of zero resonance width, the two fermions cannot penetrate the wall 
of infinite repulsion that holds them to geth er. The wave function of such a state is non-zero for only a few relative 
separations. As was shown in Section [I A , the relative motion of two particles with total Bloch momentum P is 
equivalent to the motion of a single particle on the same lattice with hopping amplitudes t'{r) = 2<(r) cos (P-r/2) in 
the external potential (?(r). As the total momentum departs from 0, the variation of the hopping terms is quadratic in 
P. The second derivative of the energy with respect to P gives the effective mass. Assuming that the wave function 
of the stable state ip{r) is non-degenerate at P = 0, we treat the variation in the kinetic energy for small P as a 
perturbation. To lowest order. 



E{P) = £:(0) - ^ - rOV-W 



p. (r _ r') 



(5.11) 



from which one obtains the mass. To this order, the mass is finite only for states with ip{r) 7^ on more than one 
site. When the state becomes unstable, its energy is altered. Since A(r) is an implicit function of the total Bloch 
momentum (through the variation of the fermion kinetic energy), Equation (^.5| ) contains the next-order correction 
to the energy. This correction is proportional to the inverse height of the repulsion walls and is therefore larger than 
the linewidth ( 5.10 ). Equation (|5.5| ) determines the bosonic mass when ( 5.11 ) does not disperse with P, as in the 
case of resonances formed on one single site. 



B. Resonance on a square lattice: an example 



Here we illustrate the results from the previous section. Consider a resonance that is formed on a square lattice by 
two fermions on the same site. This is possible when the hopping range is restricted to nearest neighbors (amplitude 
t) and the inter-particle potential is large for two fermions on adjacent sites, gi ^ \t\. The value of the potential go on 
the same site controls the energy of the resonance and therefore should be within the continuum of two free fermions: 
-8\t\ <go<8\t\. 

We start with a stable state, which is formed in the limit of gi 00. In this case, the bound state has the wave 
function ip{r) = 5ro and the energy E = go. The value of A(r) = g{r)ip{r) on sites with an infinite repulsion is then 
calculated by standard perturbation methods (with ( |2.7| ) as a perturbation Hamiltonian) : 

(5.12) 



A(r) = lim ,gi?A(r) = -2t cos (P-r/2). 



We then have a 5 x 1 column vector A, 



A = 



/ A(0,0) \ 
A(1,0) 
A(0,1) 
A(-1,0) 

VA(o,-i)y 



9o \ 

-2tcos {Pxa/2) 
-2t cos {Pya/2) 
-2t cos {Rj; a/ 2) 
-2tcos{Pya/2) J 



(5.13) 



For spin-singlet states, we only need that part which is symmetric under inversion of the relative separation r. 
Performing a unitary transformation, we obtain a 3 x 1 column vector for singlet states. 



A 



90 "o(P) 
-2t lii(P) 

-2t U2(P) 



(5.14) 
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where 



uo{P) = 1, 
ui(P) = cos(P^a/2) + cos(Pya/2), 
U2(P) = cos {Pxa/2) - cos (Pya/2). 



(5.15) 



In the hmit as gi — > oo, the matrix g ^ is given by 
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9o 





(5.16) 



A large but finite repulsion for nearest neighbors gives 















91' 











9^' 



(5.17) 



We then obtain the following expression for the energy of a resonance through the first order in t/gi: 



= 90 



Af2 

— K(P/2)+^2(p/2) 

gi 



This gives the inverse effective mass of the resonance boson as 

1 a2 



and the width of the bosonic band as 



AE- 



gi 



gi 



To determine the resonance width, we use Equations. ( 5.101 ) and ( |5.9| ) with 

2^ /2t[ul{P/2) + ul{P/2)] 
9oUi{P/2) 
5o"2(P/2) 



gogi 



(5.18) 



(5.19) 



(5.20) 



(5.21) 



Evaluation of d(j)^Ad(p is com plicated by the fact that its inverse, A ^ — g ^ — G, has a zero eigenvalue. Indeed, it 
follows from Equation ( 2.22 ) that 



^-^A = 0. 



(5.22) 



However, by its definition, vector d(j) is orthogonal to A and d(f>^Ad(l) is not singular. The calculation simplifies 
considerably for zero Bloch momentum of the resonance: U2(0) = and the linear space becomes two-dimensional. In 
this case, d(j) is also an eigenvector of the matrix A~^, with eigenvalue A. The set of four equations. Equation (5.22) 
and A~^ dif) — X d(f), contains three identities for the matrix Q (see Appendix |^) and determines the eigenvalue A. 
Solving for A yields 



A = - 



go 



go 

At 



no- 



(5.23) 



Upon substituting the expression (5.23) into Equation (5.1C) and expressing Qiq in terms of the simpler quantity GoO) 
we finally obtain 



ImE = 



me 

51 



00 



00 1 



(5.24) 



where Qqq = Qm{E ± «0+,0). We see that, for both choices of the sign, the position of the pole is shifted onto the 
unphysical sheet. The resonance width is second-order in t/ gi, while the width of the boson band (5.20) is first-order. 
Therefore, in the limit of a small i/ffii resonance bosons can propagate during their lifetime. 
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VI. CONCLUSION 



We have discussed the properties of extended two-fermion bound states and resonances on a lattice by using an 
idealized Hubbard-like model with finite-range interactions. This two-body problem is exactly solvable on any periodic 
lattice, so that all approximations are well controlled. The two-body states may be relevant to high-temperature 
superconductors for two reasons. First, when the size of a Cooper pair in a superconductor is comparable to the 
average inter-particle distance, a two-fermion picture may provide a better starting point than a Fermi-sea description. 
The superconducting transition would then be similar to the Bose condensation in liquid ^He. Second, any description 
of c?-wave superconductivity must go beyond a zero-range interaction between fermions — such as in the BCS theory 
— simply because the d-wave order parameter vanishes for two fermions on the same site. 

The symmetry of the order parameter A(r,r') in our picture is related to that of the ground state of a pair. This 
observation pinpoints the important role of the fermion kinetic energy in determining the symmetry of the energy 
gap. It appears to us that the d-wave symmetry is related to the reversed sign of the hopping amplitudes along the 
diagonals of Cu~0 sheets, a consequence of the c?a;2_j,2 character of doped holes on copper sites. We have demonstrated 
that, when the magnitude of such an amplitude is comparable to that for nearest neighbors, the ground state of two 
fermions will be a dr^2_y2 wave. 

In the strong-binding case, we have shown that pairs formed by fermions on adjacent lattice sites can easily move 
through the lattice. Thus, the condensation temperature need not be vanishingly low in this limit. This observation 
resloves the large-effective-mass problem raised by Nozieres and Schmitt-Rink in this connection |Q] . Our conclusions 
in this regard are valid for narrow two-fermion resonances as well. 

For weakly bound fermions, we have found that the size of a pair having a non-trivial symmetry tends to be smaller 
than the size of an s wave. Only in the latter case is the size determined by the depth of the bound-state energy 
level and given by the tunneling length = fivp/^- In the former case, however, a pair's size is also sensitive to a 
short-range cutoff, which we have identified with the range of the attractive interaction. This property has its origin 
in the weaker tunneling exhibited by states with higher angular momenta. It is therefore not surprizing that the 
cuprates, quite possibly d-wave superconductors, have remarkably small Cooper pairs. 
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APPENDIX A: NEAREST-NEIGHBOR ATTRACTION AND MANY-BODY BOUND STATES 



To clarify the issue of many-body bound states raised in Section III, we consider the following simple model on 
a square lattice. In the limit of strong attraction the fermion kinetic energy can be neglected. The interactions are 
given by an infinite on-site repulsion, an attractive potential gs < for nearest neighbors in the spin-singlet state, 
and a repulsive potential gx > for nearest neighbors in the spin-triplet states. Then for two nearest neighbors, the 
Hamiltonian can be written as 

^^^n.n. = ^■^^^■^■^ + (^j, _ 55) (Si -82) 

_ gr + gs gr-gs ^ 

where Si and S2 are fermionic spins and P12 is the exchange operator. 
Three fermions on a square lattice will lower their energy to 

E1+1 = 5S < (A2) 

by forming a pair and leaving one fermion free. Alternatively, there may be two bonds connecting one fermion with 
the other two. The lowest energy of this system will be given by the smallest eigenvalue of the Hamiltonian, 

g3 ^ t^^^ + (gT - gs)Si-(S2 + S3), (A3) 
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which is 



Es = (A4) 



The "2+1" configuration then has a lower energy if gr > —gs > 0. This condition can be relaxed in a ferniionic 
ensemble. Indeed, if there are two sets of "2+1" clusters on the lattice, the two free fermions will pair and reduce the 
energy further. Therefore, in a large ensemble, 

E2+1 = ^, (A5) 

and the stability of pairs requires that gx > 0, i.e. the triplet interaction must be repulsive. 

Four fermions have a larger number of states from which to choose. The paired state "2+2" has the lowest energy 

E2+2 = 2gs. (A6) 

Other possibilities include clusters of four spins having three or four bonds [Figures §(a), (b), and (c)] with the 
following energies: 

) = (3 + V3)g5 + (3 - V3)gT ^ ^^^^^ 

i?f=2.95 + 5T, (A7b) 
4'^=35s + 5T. (A7c) 



A comparison between Equations (A6) and ( A7a , b, and c) shows that pairs will merge into squares unless there is a 
strong enough triplet repulsion gx. We thus obtain a condition necessary to prevent the many-body clustering in the 
limit of strong attraction: 

9t > -9s- (A8) 

Extending the analysis in this form to larger clusters of fermions is impractical due to an exponentially growing 
number of possible ground states. Instead, we shall give an upper bound for the triplet repulsion that ensures the 
stability of pairs against the formation of many-body bound states. The idea is based on the following observation. 
The ground-state energy of a fermionic cluster with N bonds is obtained by minimizing the expectation value of the 
Hamiltonian 



+ (5T-.gs)(S.-S,)j. (A9) 

Minimizing the energy of each bond separately gives a crude estimate of this energy as Ngs, obviously much lower 
than it is. One can do better, however, by minimizing the energies of separate small clusters of bonds (not spins). 
Indeed, the lowest energy of a small bond cluster can also be found by minimization over all states of the large spin 
system. It is clear then that such a separate minimization results in a lower energy. If the bonds are grouped in k 
small clusters, we have 

H = Hi+H2 + ... + Hk (AlO) 

and therefore 

min(iJ) > min(i/i) + min(iJ2) + . • . + min(i/fe). (All) 

The equality holds if all minima are reached in the same state of the spin system. In general, however, [Hi,Hj\ ^ 
(whenever clusters i and j have common spins) and thus the eigenstates of Hi and Hj do not coincide. This 
consideration allows one to obtain a reasonable lower bound for the ground-state energy of any cluster having an even 
number of bonds. It can be shown that the Hamiltonian is a sum over pairs of bonds (with the two bonds of any pair 



sharing a common spin). The ground-state energy of two such bonds is E3 given above Equation (A4)]. Therefore, 
the energy of a cluster with Nb bonds satisfies the inequality 

En. > m'-^^. (A12) 
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In the dimerized state, Ns/2 spin pairs will have the energy 



En, = 



9s_ 
2 ' 



(A13) 



which should be lower than _BjVi, to prevent the many-body formation. Replacing _EjVi, by its lower bound (A12), we 
find that 



gr > -gs 3 - 



2N, 



(A14) 



is sufficient for the stability of pairs. The lowest possible ratio for Ng/Nt, on a square lattice is 1/2, achieved when all 
sites on a large lattice are occupied by spins. Therefore, 



gr > -2gs. 



(A15) 



provides a sufficient condition for a square lattice. 

One can see, however, that this condition can be relaxed. As has been mentioned, the largest number of bonds 
originating from one spin is achieved when all spins have four nearest neighbors. In this situation, one improves the 
lower bound for the energy by dividing the system into clusters with four bonds originating from one spin [Figure 
^(d)]. Thus, (A15) can be replaced by a milder condition. 

A cluster with an odd number of bonds can be divided into a cluster with an even number of bonds and another 
one with three bonds. Then one can readily verify that condition ( A15| ) again ensures the stability of pairs. 



APPENDIX B: POINT GROUPS OF SQUARE AND CUBIC LATTICES 

The point group of a square lattice consists of eight elements: a unit element; three proper rotations about the center 
of the square (90, 180, 270 degrees); two reflections with respect to the coordinate axes that are parallel to the square's 
sides and that have the square's center as the origin; two reflections with respect to the two diagonals of the square. 
These eight elements are isomorphic to the group of matrices ±/2 , ±ti , ±zr2 and ira with I2 the 2x2 unit matrix and 
Tj the Pauli matrices. There are five conjugate classes, which implies the following five irreducible representations: 
a unit representation Ai; three non-trivial one-dimensional representations A2, Bi, and B2, implemented by the 
functions xiy2 — X2yi, — , and xy, respectively; a two-dimensional representation E, transforming as the pair of 
fimctions (x, y). 

The proper cubic group consists of 24 elements. In addition to the unit element, there are three C4 axes running 
through the centers of opposite faces of the cube (9 rotations); four C3 axes along the diagonals (8 rotations); and 
6 C2 axes connecting opposite edges (6 rotations). These elements form five conjugate classes, which implies the 
existance of five irreducible representations: two one-dimensional representations Ai and implemented by the 
functions 1 and xyz; a two-dimensional representation E transforming as the doublet (x^ — 2z^ — x"^ — y"^)] and 
two three-dimensional representations Ti and T2, implemented by the triplets {x,y,z) and {yz, zx,xy). 

APPENDIX C: THE GREEN'S FUNCTION AT P = FOR A SQUARE LATTICE 



The two-particle Green's function was introduced in Section [IB. At a total Bloch momentum P = 



GiE; r, r') . G^E, 0; r, r') ^ 1 ^ (CI) 



E~E„ 
q ^ 



where ■(/'(r|q) = a/2 cos (q-r) for r 7^ 0, V'(0|q) = 1; and -Eq is the kinetic energy of two fermions with Bloch momenta 
q and — q, 

Eq = 2j2ti^y^"'- (C2) 

r 

Introduce 
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N ^ E-E„ 
q ^ 



which is the Green's function for unsymmctrizcd states, with '0(r|q) — e^"^'^ . If r 7^ and r' 7^ 0, 

G{E- r, r') = ij{E- r - r') + ipiE; r + r'), (C4) 
G(£;;r,0) = G(i;;0,r) = V2V'(-B;r), (C5) 
G{E:Q,Q)^tl){E-Q). (C6) 

It is straightforward to verify that ?/'(£'; r) satisfies the equation 

Ei:{E- r) - 2 ^ t{v')i^{E- r + r') - <5r,o. (C7) 

r' 

We shall restrict the hopping range to nearest neighbors (amplitude ti) and next-nearest neighbors (amplitude ^2)- 
Then the kinetic energy of two fermions is 

Eci ~ Ati{cos q^a + cosgj,a) + 8t2 cosqxacosqya. (C8) 



The minimum of iJq appears at q = when both ti < and ti + 2^2 < 0. As has been shown in section II C, the 
ground-state wave function of two fermions has an alternating sign only if there exist positive hopping amplitudes. 
We therefore consider the case <2 > 0, which may be relevant to cuprate superconductors. 

At P = 0, the two-fermion energy band Hes between i?min — -E'(o,o) — 8^1 4- 8t2 and £^max = -E(Tr,7r) = — 8ti + 8^2- 
The density of states diverges when the line of constant E touches the boundary of the Brillouin zone, at i?vH — 
^(tt.o) = —8^2- It will be convenient to measure the energy from the van Hove singularity: 

E' = E + ?>t2, (C9) 
ELn = 8^1 + 16^2 < 0, (CIO) 
E[^^^^ = -8ii + 16^2 > 0. (Cll) 

On an infinite lattice, the Green's function for relative separation r = (m, n) is 

i/'(m,n)= (C12) 
d^q cos mqx cos nqy 



(27r)2 E — Ati{cosqx + cosqy) — 8^2 cosqx cosqy ' 
When n = 0, the integration over qy is particularly simple. We have 



V'(m,0)== (C13) 
1 cosmqx dqx 



[E' - E'^^,^^ cos2 {qx/2)] [E' - E'^^^ sin' (fe/2)] 

The sign of the square root is determined as follows. When the root is real, it has the same sign as £", which ensures 
that V'(^';0,0) 1/E as E ^ 00. The region for imaginary values of the square root, between the branch points 
E'^^^sin' {qx/2) and E^^^cos'^ {qx/2), is accessed either from E' = E'^-^^^cos'^ {qx/2) + 0+ > (by encircling the 
singularity counterclockwise), or from E' = E'^^^sin' {Qx/'^) — 0+ < (by encircling the singularity clockwise). In 
both cases, the square root has a positive imaginary part and zero real part. 

At the origin, ijj can be expressed in terms of a complete elliptic integral of the first kind, 

^(0,0) = ^K(fc), (CM) 

where 

fc2 = 1 - (1 - E'^,jE')il - E[^JE'). (C15) 
This function can be analytically continued into the region ii'jnin < E' < E'^^^. using the identity 
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^ K(iA)±.K-(iA) 

k 

(the sign is chosen in accordance with the rule stated above) . The density of states is then given by 

-'-)-^^^'<'<'/''' (-) 

This function is plotted in Figure ^ for ti = —1 and several positive values of t2. As t2 —ti/2, the band flattens 
at the bottom and the logarithmic van Hove singularity coalesces into the branch point of the Green's function. This 
strongly enhances the density of states near the lower edge of the band. 

As m increases, expressions for ipiE', m, 0) become progressively more complicated and include, in general, complete 
elliptic integrals of all three kinds. For example, 

where the integral of the third kind is defined as 

U{a\k)^ C ^ '^'^ =. (C19) 



An expression like ( |C18| ) provides little insight into the behavior of the Green's function. Fortunately, it is possible 
to obtain an approximation for ^{E; r) in the vicinity of its singular points. Several authors [p^ have identified 
the enhanced density of states near the van Hove singularity as a possible candidate for explaining high critical 
temperatures in cuprate superconductors. In view of that, we discuss the leading behavior of ^{E; r) near E — E^h, 
i.e. E' = 0. 

As E' ^ 0, k ^ oo, so that K(l/fc) 7r/2 and K'(l/fc) - ln4fc ^ 0. Therefore, near E' = 0, 

Note that Re-!/'(0,0) ±1/Eq as E' ^ 0^, in accordance with the sign rule for the square root. At even sites along 
the axes, the difference ip{2m, 0) — ?/'(0, 0) remains finite as E' 0: 

i/'(2m,0) - i/'(0,0) — 7 / dq^ 



ttEIj Jo sinqj; 

= (C21) 



In a similar way, 

V'(2m + 1,0) - V^(1,0) ^ 0. (C22) 

It is also seen from ( C13| ) that the two divergent contributions to 'ip{2m + 1, 0), near = and q^ = tt, have opposite 
signs and thus tend to cancel each other. In the limit as E' — > 0, 

^(i'0)-^4-r5?7iog4^- (C23) 

Eq ttEq --fcmin 

This shows that tp{2m + 1,0) is not singular at i? = i?vH- 

Using Equation (C7) it is now possible to obtain approximations for tp{r) on all sites. We demonstrate how this is 
done for the singular part of the Green's function. As i? ^ E^n = — 8t2, Equation (C7) reads 

- 8i2^(r) - 2ti J2 V'(r') - 2h ^ ^(r") = 0{1), (C24) 
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where r' refers to the four nearest neighbors of r, while r" refers to the four next nearest sites. Starting with r = 0, 
we obtain (by using symmetry of the lattice) 

^'(la) = -VX0,0) + C'(1). (C25) 



Applying Equation (C24) to site (1,0), we find that 

V'(2,l) = 0(l), (C26) 

because the divergent contributions of sites (0,0) and (2,0) are cancelled by those of (1,1) and (1,-1). We now 
separate the lattice into two sublattices, putting adjacent sites on different sublattices. In general, ip{r) is singular on 
the sublattice that contains the origin: at sites with two even coordinates, ip{r) = -0(0,0) + 0(1); at sites with two 
odd coordinates, ^(r) = —0(0, 0) + 0(1). On the other sublattice, ip{r) is regular at £' = E^n. This means that the 
relative motion of the particles at E = E^h is restricted to one of the sublattices. 



As we mentioned, on any site r, "0(1") is a linear combination of three complete elliptic integral with modulus ( C15 ) 



Therefore, it should be possible to find "0(1") on the whole lattice if it is known on three conveniently chos en s ites. For 



example, one starts out with (0,0), (1,0), and (2,0). Then -0(1,1) is found from the Laplace equation (C7) around 
site (0,0) by using the symmetry arguments ^{m,n) = ip(n,m) and il){—m,n) = ^j{m,n). The Laplace equation on 
site (1,0) expresses ?/'(2, 1) in terms of -0 at sites (0,0), (1,0), (1,1), and (2,0). To proceed, we use an idea due to 
Friedberg and Martin Q, who showed that there is a lattice analog of the gradient condition, VV'(r)|| r. Extending 
their idea to the case of several hopping amplitudes, we find an additional relation between values of '0(r + r') on sites 
surrounding r. 

5:'^^(r + r') = E'^^(r + r'), (C27) 

r' r' 

where means that sites with x' = {y' = 0) are excluded. On sites with high symmetry, such as {m,m), (m, 0), 
and (0,m), this equation becomes an identity and supplies no additional information. 

Now we are in a position to recover "0 on the remainder of the lattice. We set up the Laplace equations on sites 
(2,0) and (2,1) in order to relate at sites (3,0), (3,1), and (3,2) to each other and to V' at sites where it has been 
calculated. The gradient equation on site (2,1) becomes a third one in this set of linear equations, which yields three 
more sites. V'(3, 3) is obtained from the Laplace equation on site (2,2). This establishes the algorithm for computing 
ip on the remaining sites. When ip is known on the square with the corners (0,0), (0, m), (m,0), and (m,m)^ we 
write the Laplace equation for site (m, 0) and the Laplace and gradient equations for site (m, 1). We then solve these 
equations for at {m + 1, 0), (m + 1, 1), and (m + 1, 2). Using the Laplace equations for sites (m, 2), . . . , (m, m), we 
move up and complete the side of the square. 

In two particular cases, ^2 = and 2^2 — —ti, the calculation of 0(r) is simplified by the absence of elliptic integrals 
of the third kind. In the latter case, the van Hove singularity coalesces with the bottom of the continuum band, 
^min — resulting in a stronger divergence of ip{E] r) near E' — 0. We write down an approximate expression for 
'ip{E,r) without presenting a derivation. Just outside the continuum band, as E' ^ 0^ , 

1 (C28) 

16K' 



^\ E' / ^ 2fc - 1 ^ 2/ - 1 



fe=i 1=1 

which can be readily continued to the region Q < E' < E'^^^; the point E' ^ should be encircled clockwise. 

APPENDIX D: PROOF OF IDENTITIES USED IN SECTION V 

To prove the relation 

-a.|a.i. (D1) 

we write the quantity —A^{dG/dE)A = d in its explicit form, 
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(D2a) 



hj=0 



It follows then from the definition of the Green's function (2.21) and the orthogonality relation (2.14) for functions 
ij{r\p) that 



The sum must be taken over the entire latti ce, s o that this identity cannot be applied directly to the nx n matrix 
0{E,P). We substitute this expression into ( D2a ): 



'^^ E E^*(I^0ff(Ri)G(Ri,r)G(r,R,)5(R,)V'(R,) 

i,j—0 r 

= V'*(r")5(r")G(r",r)G(r,r')g(r')V(r') 

r,r',r" 

= ^r(r)^(r) = 1. 



(D2b) 
(D2c) 



In ( D2b ), the n ewly add ed te rms do not contribute bec ause either V'(r^ ) = or ^(r') = for sites r' ^ Hj. The 
transition from ( D2b ) to ( D2c ) is justified by Equations ( 5.1a ) and ( 5.1b ). 

Next, we prove that, in A{dQ /dE)A, DQ/dE can be replaced by —Q^. Explicitly, 



X: A(R.,R,)^^^^^5(R.)^(R.) 



].k=0 



J2 ^A(R.,R,)G(R,,r)G(r,R;.)3(Rfe)^/.(Rfc) 

n 

-^5]A(R„R,)G(R„r)^(r) 

3=0 r 



(D3) 
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(a) 



(b) 



(c) 



FIG. 1. Deep bound states on an f.c.c. lattice. Symmetries and energies of the states: (a) Ai, gi + 8t. (b) E, gi — At. (c) 
T2, Si. 
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FIG. 2. Shallow s-wave and d-wave ground states on a square lattice. Nearest-neighbor hopping = — 1. 
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FIG. 3. Bond clusters on a square lattice. 
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FIG. 4. Density of states for different values of nearest-neighbor hopping t2. ti = —1. Values of t2 are shown on the graph. 
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